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SPECTRUM OF TURBULENCE IN A CONTRACTING STREAM 1 

By H. S. Ribner and M. Tucker 


SUMMARY 

The spectrum concept is employed to study the selective effect 
of a stream contraction on the longitudinal and lateral turbulent 
velocity fluctuations of the stream. By a consideration of the 
effect of the stream contraction on a single plane sinusoidal 
disturbance wave, mathematically not dissimilar to a triply 
periodic disturbance treated by Q. I. Taylor, the effect on the 
spectrum tensor of the turbulence and hence on the correlation 
tensor is determined. Lack of interference between waves 
; follows from the postulation of a low level of turbulence; this 
and the assumption of an invisdd fluid imply neglect of decay 
effects. The compressibility of the main stream is taken into 
account, but the density fluctuations associated with the turbulence 
are assumed to be negligible; this would be the case if the turbu- 
lence originated from wakes and boundary layers in the very 
low-speed portion of the flow. For an axisijrnrnetric contraction 
and a particular isotropic initial turbulence, some explicit 
results are obtained. The one-dimensional longitudinal 
spectrum is found to be. distorted {as well as reduced in ampli- 
tude) with its peak shifted well to the right of the initial position 
above the zero of the wave-number scale. The selective effect of the 
contraction on the mean square longitudinal and lateral com- 
ponents of turbulent velocity is found to be given uniquely when 
the initial turbulence is isotropic, regardless of the details of 
the spectrum. If the initial spectrum is anisotropic, as, for 
instance, that produced by a damping screen, then the selective 
effect is altered. 

In a crude extension, decay effects outside the scope of the 
theory are allowed for in first approximation. With this 
extension, a comparison with experiment is made of the selective 
effect on turbulent intensity where the estimated decay effects 
are comparable with the contraction effects. The agreement is 
good for the longitudinal component, very poor for the lateral 
component, the experimental data themselves being in conflict. 

INTRODUCTION 

The generation of a wind-tunnel flow is always accom- 
panied by a certain amount of turbulence; this is one respect 
in which the flow fails to simulate free-flight conditions. 
Measurements in the tunnel, particularly those sensitive to 
boundary-layer behavior, are known to be affected by this 
turbulence. Accordingly, the tunnel designer attempts to 
reduce the intensity to the lowest practicable level. The 

' Supersedes NAOA TN 2006, "Spectrum of Turbulence In a Contracting Stream," by 
H, S. Rlbner and M. Tucker, 1062. 


use of honeycombs and damping screens in a large low-speed 
section (settling chamber) followed by a sharp contraction 
to the much-higher-speed working section is known to be 
effective. The honeycombs and screens located in a low- 
speed section reduce the absolute level of the turbulence 
with little drag penalty; then the relative level is greatly 
reduced by the large gain in tunnel speed through the con- 
traction, aside from any effect of the contraction on the 
absolute level. 

Once the characteristics of honeycombs and screens are 
known, the further quantitative estimate of the reduction in 
turbulence involves a knowledge of the effect of the tunnel 
contraction 2 on the turbulence. It is known that the 
longitudinal component of the turbulence is greatly reduced 
(in absolute value) by the contraction; the behavior of the 
lateral component appears, on the other hand, to vary from 
no change to a substantial increase. Prandtl (ref. 1) ob- 
tained a quantitative estimate of the first effect by consider- 
ing the conservation of energy for a perturbed longitudinal' 
filament: if the initial stream speed is U, the filament speed 
TJ-\-u, with tt-C TJ, and the final stream speed is IU, then 
the final filament speed must be lU-\-l~ l u; that is, the con- 
traction reduces the longitudinal perturbation velocity -a, by 
the factor Z _1 . For the lateral effect, Prandtl applied con- 
servation of momentum to a small rotating cylinder of the 
fluid, with its axis cross stream, as the fluid traversed 
the tunnel contraction. He concluded that the lateral 
perturbation velocity v is increased by a factor fl. 

Prandtl’s considerations on the effect of a stream con- 
traction were limited, as has been noted, to particular 
idealized “eddies.” G. I. Taylor (ref. 2) later attempted 
more realism by treating a mathematically defined model of 
turbulence which amounted to vortices in parallelepiped 
partitions arranged in a regular three-dimensional array. 
The changes in vorticity on traversing the contraction were 
determined from a theorem based on conservation of circu- 
lation for an inviscid fluid ; the corresponding altered turbu- 
lent yelocity pattern was then calculated. The final result 
of the analysis consisted in expressions for the root-mean- 
square longitudinal and lateral turbulent velocity compo- 
nents u' and t/ downstream of the contraction expressed as 
ratios of the corresponding values upstream. 

• 

5 The considerations of this paper are not limited to a wind-tmmel contraction: they may 
be applied to any stream tube of varying cross section largo compared with the scale of the 
turbulence. 
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The initial condition of isotropic turbulence (mean values 
unaffected by rotation or reflection of axes) was approxi- 
mated by specifying the vortex partitions to be cubical. For 
this case the reduction in the longitudinal component u' was 
found to vary more nearly like 1.5 l~ l than the value Z -1 sug- 
gested by Prandtl. No explicit result was found for the 
variation of the lateral component, however: the calculations 
contained a free parameter. 

Taylor’s results for the longitudinal component agreed 
fairly well with the experimental data then available, but it 
is now considered that the measurements were made too dose 
behind the screens for the screen-produced turbulence to 
have been isotropic. On theoretical grounds, the objection 
to Taylor’s theory is threefold: first, the decay processes of 
turbulent mixing and viscous dissipation, which result in a 
reduction of the mean intensity with axial distance in the 
wind tunnel, are neglected; second, the assumed model of 
turbulence fails to exhibit the spatial and temporal random- 
ness of actual turbulence; third, no choice of the parameters 
in Taylor’s model corresponds to isotropy. In a sense all 
three objections apply likewise to Prandtl’s results: no model 
was employed in his considerations, and hence no distinctions 
between the effects of isotropy and anisotropy were made. 

The second and third objections can be removed by work- 
ing, not with a model of turbulence, but instead with a 
Fourier integral representation of a random turbulent field. 
The integral can be interpreted as a superposition of plane 
transverse sinusoidal waves of all wave lengths and with 
apparently random orientations and phases. This aggregate 
of plane waves constitutes the (three-dimensional) spectrum 
of the turbulence. Only the statistical aspects of this 
spectrum will be known, not, for example, the detailed phase 
relations. Mean-square velocity components 'may be ob- 
tained by an integral of certain spectrum functions in which 
the phase relations are suppressed; these functions are in- 
cluded in the ‘spectrum tensor’ (ref. 3). 

Taylor’s concepts may be applied to find the effect of a 
stream contraction on a single plane wave. The effect 
under the assumptions is linear; therefore, the superposition 
implied by the Fourier integral may be employed to obtain 
the contraction effect on a field of turbulence. In particular, 
if the initial spectrum tensor is known, the final spectrum ten- 
sor is determined; the initial and final mean-square velocity 
components then result from quadratures. Indeed, from 
the same information, changes due to the stream contraction 
in correlations of velocity at different points may be calcu- 
lated: use is made of the fact that the correlation tensor 
is an inverse Fourier transform of the spectrum tensor 
(ref. 3). 

Accordingly, the injection of the spectrum point of view 
into Taylor’s original concept of the contraction effect 
makes possible a more realistic calculation of the changes 
in mean-square velocity components. In addition, it 
provides much more detailed information concerning 
changes in the statistical properties of the turbulence; that 


is, it provides the changes in the spectrum tensor and in the 
correlation tensor. 

The ideas just outlined are developed in the present paper. 
The first section is devoted to an account of turbulent 
spectrum analysis in a form specially adapted to the analysis 
of the contraction effect. In this account, which is a general- 
ization of a development in reference 4, the role of the spec- 
trum tensor is subordinated to that of the individual Fourier 
components (plane waves) in contradistinction to the cus- 
tomary treatment. This approach has perhaps an auxiliary 
merit in providing some better physical insight into the 
significance of the spectrum tensor. 

Next the effect of a stream contraction on a single plane 
wave is calculated by an applicatipn of Taylor’s concepts. 
The treatment is slightly more general in that compressibility 
of the main stream is allowed for. The density fluctuations 
associated with the turbulence are assumed to be negligible; 
this would be the case if the turbulence originated entirely 
from boundary layers and wakes in the very-low-speed portion 
of the flow. Following Taylor, the problem is linearized by 
postulating a sufficiently weak turbulence so that the self- 
distortion of the turbulent eddies is small compared ■with the 
distortion imposed by the contraction of the main stream; 
this together with the assumption of an inviscid fluid implies 
neglect of the decay of the turbulence. 

In succeeding sections the spectrum and correlation tensors 
downstream of the contraction are expressed in terms of 
the corresponding initial tensors. For the special case of 
an axisymmetric contraction and isotropic initial turbulence, 
the ratios of the root-mean-squaro longitudinal and lateral 
Velocity fluctuations downstream and upstream are obtained 
explicitly in terms of the parameters defining the contraction. 
For a particular subcase where the initial isotropic spectrum 
tensor is specified, the corresponding ‘one-dimensional’ 
spectrums (as would be recorded by stationary liot-wiro 
probes) upstream and downstream of the contraction are 
calculated; the specification is such that tho upstream one- 
dimensional spectrum corresponds to experiment (ref. 6, 
p. 35) in a number of cases of isotropic turbulence. 

Most of the calculated contraction effects are amenable to 
experimental checks either directly or indirectly. The 
available experimental data, however, are limited to tho 
changes in the root-mean-square velocities. A comparison 
with these experimental data is given with an estimated 
allowance for decay effects' outside the scope of the theory. 
Design curves of the changes in the root-mean-square velocity 
components neglecting decay are included for engineering 
purposes. 

SPECTRUM ANALYSIS 

Representation of turbulence by superposition of plane 
sinusoidal waves. — Suppose q h q it g s represent tho com- 
ponents of velocity in a turbulent field; that is q_ h g a , and g 3 
vary in an apparently random manner in space and time, 
and the mean values 2 i= 2 s= 2 j= 0 . Subject to certain 
conditions, a snapshot of this field at any instant can bo 
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represented as a set of three-dimensional Fourier integrals 

co . 

<la (xi, Xu, x *) = fffQ a (*«> ki ’ * 3 ) e Hk ^+ k ^ k ^> dhdktdkj 

— CD 

where a=l, 2, or 3 and the significance of k h k 3 , and k 3 
will be brought out later. A continuous representation of 
the turbulent field is obtained by allowing the Q a to vary 
with time. 

It will be convenient to abbreviate the Fourier integral to 

CD 

ia(x)= j J j Q a (k) e {t -*-dr (k) (la) ■ 

— co 

and to introduce the companion equation 

Q a (k)=8^J J jq a (x)e- { ^r(x) (lb) 

— CD 

whore ]c=k 1 , k 2 , k 3 ; x=x } , Xi, x 3 ; dr (x)=dx l dx 2 dx 3 . The 
second equation allows, in principle at least, the coefficients 
QAb) to be calculated. Mathematically, g a (x) and Q a (k) are 
termed three-dimensional Fourier transforms of each other; 
use will be made of this relation later. 

The velocity components q a are connected by the condition 
of continuity. In many cases of practical interest, these 
turbulent velocities originate from boundary layers and the 
wakes of obstacles in flows of low subsonic speed, so that 
associated density fluctuations may be ignored; this is still 
permissible when the turbulence so produced is transported 
by a high-speed stream. Thus the incompressible form of 
the continuity equation may be used and the result is 

Qiki -{- Qtk-i -f - Q 3 k 3 = 0 

This relation may be written more compactly as 

'22Qaka=0 ( 2 ) 

a 

Physical interpretation. — The amplitude components Q a 
are complex in general. According to equation (lb), then, 
the requirement that the velocity components g a be real 
implies that Q. a (—k) is the complex conjugate of Q,a(k). If 
corresponding terms for k and —k in equation (la) are paired, 
then sum is thus equal to the real quantity 

2(Re Q a ) cos (Jc-x)— 2(/m Q a ) sin (k-x) (3) 

The imaginary parts cancel in the pairing, which implies that 
they contribute nothing to the integral. Expression (3) 
represents a pair of plane standing waves, a cosine wave and a 
sine wave, with normals in the direction k=(k lt k 2 , k 3 ), 
where x=(xi, x 2 , x 3 ) is the radius vector to any point. The 
vector k is termed the wave-number vector and its magnitude 
k simply the wave number; the corresponding wave length is 
27r divided by the wave number. Since k is perpendicular to 
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the wave front, it is sometimes referred to herein as the Svave 
normal.’ 

The continuity condition, equation (2), states that both 
the real part (Re Q a ) and the imaginary part (Jm Q a ) of the 
amplitude vector Q=Q a = (Qi, Qi, QA are perpendicular to 
the wave nbrmal, k ; that is, both waves of expression (3) are 
transverse. For each wave any one of the parallel planes 
containing both the local velocity vector q and the wave- 
number vector k (which is perpendicular to Q and hence to_g) 
is called the plane of polarization. The cosine wave (real 
part) and sine wave (imaginary part) may be polarized in 
different planes in general; the necessary and sufficient con- 
dition that they be polarized in the same plane is 

Re & _ Re Q 2 _ Re Q 3 , ... 

ImQx ImQi ImQ 3 ^ ' 

Equations (1) are now seen to represent a superposition of 
plane sinusoidal waves (Fourier components) with all orienta- 
tions of the wave fronts (all directions of the wave normal k) 
and all wave lengths (all wave numbers k). Each wave is 
transverse, and all planes of polarization, are permitted. 
For each value of k there exist a cosine wave and a sine 
wave; their respective amplitudes and planes of polarization 
are different in general. The complex amplitude components 
Qa(k) express, in their real and imaginary parts, how the 
respective amplitudes and planes of polarization vary with 
the wave-front orientation and the wave number. 

Mean values : correlation tensor. — Consider the spatial 3 
mean value of the product of the velocity component q a at x 
and the velocity component g/ at x'=x+r as x varies but 
the separation r of the two points remains fixed during the 
averaging process; this mean value is called a velocity cor- 
relation and is given the symbol Rc#(r). There are nine such 
correlations, corresponding to a=l, 2, 3, and j3=l, 2, 3. 
The form R^(r) has been shown to transform like a second- 
order tensor and has been designated (sometimes divided by 
g*) as the ‘correlation tensor’ (ref. 6). 

Evaluation of correlation tensor in terms of spectral 
quantities. — The basis of the idea that the correlation tensor 
might be expressed in terms of individual-wave parameters 
is drawn from reference 4. In that paper mean-square values 
of < 7 i, g 2 , g 3 in terms of such parameter’s are discussed; these 
constitute the diagonal terms of f? oj9 (r) evaluated at r=0. 
The following derivation of R a p(r) amounts to a generaliza- 
tion of the derivation of R n ( 0) given in that reference. 

Assume that the turbulence is confined to a large, but not 
infinite, parallelepiped of edges 2 D u 2 D 2 , 2 D 3 and vanishes 
everywhere outside. 1 The space average J?^(r) is to bo 

* If the statistical properties of the turbulence are independent of position (homogeneous 
turbulence) and time-independent, an average at a given time over all space equals an average 
at a given point (or pair of correlated points) ova* all time; a proof is given in appendix B. 
If the statistical properties vary slowly with time, the space average will still approximate a 
time avenge over an Interval just long enough to smooth out the fluctuations. 

* The space integral of the square of any component velocity is thereby bounded, which is a 
prerequisite for the existence of a Fourier Integral, equations (1); In other words, equation 
(lb) shows that Q a must depend on the volume SDiDtlh of the region within which q a differs 
from zero: if this volume Is infinite Qa Is infinite, and the Fourier integral, equation (la), does, 
not exist. 
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token for all points x in the interior, the separation r remaining constant and very small compared with D h2 , 3 - The average is 

Bjr) s jj/ jj/ 

\ — CO CO / 

or, interchanging the order of integration, 

-d^dJ'H J jf <»b<w ^ ^;? D - 8i ° ( ^y g ~ *«»<*> 

CO CO 

Put K 2 =k 3 +k 3 ', and K 3 =k z -\~k 3 . Then recalling that k=(k u k 2 , k 3 ), dr(k) =dk l dk 3 dk z , etc., and k'-r= 

kiri+h'r 3 +k 3 'r 3 


R ^ r)= rtimf J J Q&iMJ \ f ff (h{Ki-KK 3 -k 3 ,K 3 -k 3 ) 


sm_gP s sm_gP 3 dKldrKidK% j d kidk 3 dk 3 

The inner integral approaches a limit as the parallelepiped edges 2 D u 2 D it 2D 3 become indefinitely latge. Fot this 
limiting case the correlation tonsor is then 

CD 00 

^to=^ini j J J j)j) 2 2) 3 QaiKKkJQei—h.—k^—kzy-^ldkidkidk^hm f f f~ Qai&Qtf&e^-'-drik) (6) 

■ — CO — CO 

where £)] is the complex conjugate of On(k) and r is the volume 8DJ)JD 3 of the parallelepiped. 


8tt^ 

Identification of lim Q a Q* with spectrum tensor. — 


If the field of turbulence is homogeneous, the correlations 
should be independent of the volume r averaged over when this 

8-7T 3 

volume is sufficiently large; thus lim* Q a Qp* should exist. 

T »CO T 

From the mathematical standpoint, equation (5) shows that 

8tt 3 

R a p(r) is the Fourier transform of lim — - Q a (k)Qg*(k),fmd 

T— M» T 


conversely; since it is known that the integral of |F a p(r)| 
over all r is finite,^ the inverse Fourier transform relation 

ensures that lim — QaQp* exists. Finally, since the Fourier 

r i cd T 

8x 3 

transform of h function is unique, lim Qa(k)Qp*(k) may 

T — ► <» ' T 


be identified with the form Y a p(k) defined by Batchelor (ref. 3) 
as the Fourier transform of R a p(£)- The form T a p(k) is 
known as the spectrum tensor. The Fourier transform rela- 
tions connecting the spectrum tensor and the correlation 
tensor are summarized as 


CD 


r°*(*)=3^3 J J J RadR)e^ T -dr(r) 

1 — 

(6a) 

00 

R*fi(r)=J J J Tafs^e-^ldrOc) 

(6b) 


CD 


By use of the Fourier transform relation, Batchelor demon- 
strated that is a second-order tensor and obtained a 
number of its properties. Thus, for example, r a3 is complex, 
in general, with Te a =T a p*, and the diagonal elements r„ a 
are real; also, r<^(— &) = !>„(&). It is of interest to observe 


that these same properties result immediately from the 

Sir 3 

identification of lim — QJb* with r o3 . Thus Q a Qp* is 

T — ► ao T 

complex, in general; QpQ a * equals [<?«<?/)*]*; and Q a Q a * is, 
of course, real. Furthermore, since Q a (—k) = Qa*(k), 
Qa(—k)Qff*(—k) equals Q a *(k)Qi>(k); bence r a p(— k) = 

r „.(*). 

The distinction between cases where r o(3 is real and cases 
where it is complex may be given a physical interpretation. 
The product Q a Qn*, and correspondingly r a /j, is seen to be 
real when the condition equation (4) is satisfied. This 
implies that the cosine wave and sine wave associated with 
wave number k are polarized in the same plane. The alter- 
native condition where Q a Q?*, and hence r„p, is complex 
implies polarization of corresponding cosine and sine waves 
in different planes. The velocity pattern of such a pair of 
waves is quite interesting: successive velocity vectors along 
a line in the direction of the wave normal k turn progressively 
about this line in spiral fashion; the tips of the vectors trace 
out a helical curve on a cylinder of oval cross section. 

Energy spectral density. — Each of the diagonal elements 
r„, r», and r M of the spectrum tensor r a| s may be inter- 
preted as an energy spectral density. Thus, according to 
equation (6b) 


iWoW-JJJ T u (k)dr(k) 

— CD 

f 

Therefore, the differential ^ Tndr© represents the contribu- 
tion to the kinetic energy component jj Uy per unit mass 

A 

made by waves with wave number within the range dr(k ) . 



SPECTRUM OP TURBUXENCE IN A CONTRACTING STREAM 


103 


One-dimensional spectrum. — The elements of the three- 
dimensional spectrum tensor are not directly measurable; 
they may be obtained by taking the Fourier transform of the 
measured correlation tensor. A hot-wire probe placed in 
the moving stream will, however, develop a fluctuating out- 
put voltage whose (one-dimensional) frequency spectrum 
(ref. 7) is related to a diagonal element of the three-dimen- 
sional spectrum tensor. Thus by equation (6b) the contri- 
bution to the mean-square velocity component g a s (=F aa (0)) 
from all waves with wave-number components in the direc- 
tion of the a; r axis between |&i| and |i 1 |-r|d/: I | is 


F a {k{)dki = 2 


(!! 


T aa dk 3 dk 3 )d&i 


(7) 


the factor of 2 accounting for suppression of negative values 
of ki. The function F a (ki) is the one-dimensional spectrum 
corresponding to the velocity component q a ; the values 
a=l, 2, 3 correspond respectively to the longitudinal and 
two lateral spectrums. The particular spectrum obtained 
depends on the arrangement of the hot-wire probe elements. 

EFFECT OF STREAM CONTRACTION 3 

Consider now that the turbulent velocity pattern q u q 3 , q 
is carried along by an inviscid general stream with velocity 
Ufa) in the a^-direction. Consider also that q h q 3 , q 3 are 
so small that their effect on the streamlines may he neglected 
as tho flow traverses a wind-tunnel contraction. The con- 
traction will, however, distort the shape of fluid elements. 
(See fig. 1.) The vorticity distribution will be forced to 
alter accordingly to conserve the circulation about each 
element. The net result will be an altered pattern of turbu- 
lence. Each plane wave (Fourier component) Q it Q 2 , Qse*-- 
will, in fact, be altered independently under the linearizing 
assumption to be made; the, over-all effect on q h q 3 , q 3 can . 
be obtained by tho summation expressed by equation (la). 
It thus suffices to consider the effect of the contraction on 
a single representative plane wave. 

EFFECT OF CONTRACTION ON REPRESENTATIVE PLANE WAVE 

Velocity and vorticity at upstream station. — Designate by 
A a refetence station upstream of the contraction and by B a 
reference station downstream of the contraction. (See 
fig, 1.) Let a typical Fourier component (plane wave) of 
the turbulent field q a {a=l, 2, 3) at station A be represented 
at time f=0 by 

qF— Q a A e i< -- 1 (8) 

This wave, equation (8), is supposed to be carried along by 
the main stream with Velocity U. 

The vorticity u a is obtained from the curl of equation (8) 
as 

(9) ' 

P.y 

{ 0, if any pair of subscripts are equal 
1, if afty are in cyclic order 
—1, if ai 3y are in anticyclic order 



(a) 



(b) 

(a) Tnnnel geometry. 

(b) Stream tube ccomotry. 

Figdbb 1.— Schematic representation of flow contraction parameters. 

Distortion of fluid element in passage through contrac- 
tion. — Suppose the contraction is such that The stream 
velocity U is increased by a factor l t between stations A and 
B while the breadth and the height of the tunnel are reduced 
by factors k and / 3 , respectively. (See fig. 1 (a).) In 
traveling from A to B, an initially cubical element of fluid 
of edge D will be distorted into a parallelepiped of edges 
liD, IJD, HD (see fig. 1 (b)) ; a particle in the element original- 
ly (t=0) a vector distance x from a comer particle will finally 
(t=t) be found a distance £ from the comer particle, where £ 
is related to x by 

ii=hxi\ 

£3 — ^ 2 a ' 3 / > ( 10 ) 

In this argument the modification of the streamlines due to 
the turbulent velocity fluctuations has been neglected. 
This implies that the relative displacement of two adjacent 
particles due to the . superposed turbulent motion is small 
compared with the displacement due to the tunnel contrac- 
tion. This key assumption, due to Taylor (ref. 2), linear- 
izes and vastly simplifies the problem. The limitations 
imposed by the assumption are discussed later under DECAY 
CONSIDERATI ONS. 

The velocity ratio li and the lateral and vertical contrac- 
tion ratios h and h bre related by the continuity condition 

<rliHs == 1 
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where cr is the ratio of stream densities at stations B and A; 
the density is considered uniform at each station in accord- 
ance with the initial assumption of negligible turbulent 
density fluctuations. 

Vortioity at downstream station. — The vorticity is carried 
along by the flow, the fluid elements undergoing the distor- 
tion pictured in figure 1 (b), to the approximation used. 
During the motion the strength changes in s.uch a way as to 
maintain the constancy of circulation of the fluid elements. 5 
The changes are expressed by the equations for the transport 
of vorticity in the Lagrangian form, due to Cauchy (ref. 8, 
p. 205). 



8 


The right-hand side may be expressed as k • £, where 

_fa fa fa 
~ h’U’h 

defines a new wave-number vector. 

Velocity at downstream station, general case. — Tho 
velocity distribution whose curl in the &, &, £, system is 
given by equation (lib) and which satisfies continuity is 
found to be expressible in the form 

2 (12)* 

with 



where a is the density ratio between stations B and A, and 
the derivatives bfafbxp express the effect of the fluid distor- 
tion. 6 Evaluation by means of the distortion equations (10) 
yields simply 

6> a B = <rl a w a A 

or, in expanded form 

«i B = <rli(Oi A 

<jj a B = al^f 

a>3 B —<rl^} 3 A 

These equations relating downstream and upstream vorticity 
embody the entire dynamics of the contraction effect. The 
equations are not limited to the plane sinusoidal waves dis- 
cussed earlier, but apply to any (weak) vorticity distribution 
whatsoever. 

The above derivation of the vorticity changes is substan- 
tially in the form given originally by G. I. Taylor (ref. 2) for 
the case <r=l (incompressible flow). In order to assess the 
influence of the simplifying assumptions a more general 
derivation based on the Navier-Stokes equations is given in 
the section entitled DECAY CONSIDERATIONS. 

By virtue of equation (9) as applied to equation (11a), 
the vorticity at station B is obtained explicitly as 

w a B =i<rl a Jle JjCfQSeHrs (lib) 

At 

where, it will he remembered, x is the radius vector to a 
fluid particle at time t = 0 when the fluid element is at 
station A in the moving coordinate system of figure 1 (b). 
The corresponding vector to the particle at time t—t, when 
the fluid element is at station B, is £ in that figure. When 
equations (10) are used to express x in terms of £, the ex- 
ponential fax becomes, in expanded form, 

t, &i£i i fata i faZs 

-- _ ir + T" + ir 

3 This statement la exact for the postdated invlsdd fluid. The modification produced by 
the diffusive effect of viscosity, in the case of a gas, becomes appreciable for the smaller eddies 
or higher wave numbers; for this analysis a criterion for neglect of viscous effects Is 
*fc *<<| d UJdj\. (See decay considerations, equation (43).) 

• These equations refer to axes moving with some fluid particle rather than axes fixed as in 
reference 8; the form of the equations b unaffected. 



where k is the magnitude of the wave-number vector jj. 
This result is the general solution for the contraction effect 
on a single plane wave. 

Equations (12) and (13) admit of a simple (but not obvious) 
geometrical interpretation: traversal of the stream contraction 
alters the initial -plane wave, equation ( 8 ), so that -its wave- 

k\ fa k 

number vector k=(fa, fa, k 3 ) is transformed into -r-i -p 

and its amplitude vector (Q A , Q A , $/) is transformed into the 
projection of ( Q A /li , Qi A /ls, $j A /4) on a plane normal to the 
new wave-number vector k. 

Velocity at downstream station, axisymmetrio contrac- 
tion. — In case the stream contraction is axisymmetric, 7 a con- 
siderable simplification results. The condition for nxisym- 
metry h=h, with use of the continuity equation (2), reduces 
equations (13) to 


A , Qi A fa*+fa*+fa* 

” h efa’+fa'+fa* 


«•*=£ [« 


Qi A fak 3 (1 — t). 


tki iJ rfa 1J rk 3 


a , Q±fakf 1 — e) 
3 ^efa'+fa'+fa* 


(14) 


where e =lt/lf. The considerably greater complexity of 
equation (13) is perhaps obscured by the purposely expanded 
form of equations (14). 

If the initial wave normal k is perpendicular to the (longi- 
tudinal) x^axis, the component k x vanishes and equations 
(14) reduce further to 


Q B =Qi A lh' 
Qi B =Q> A lh - 
Qz B =Qffa, 


(16) 


The same equations result when Q A may be neglected in 
comparison with Q 2 A and Q 3 A , that is, when the amplitude 
vector is substantially normal to the Zi-axis. Equations (16) 
state that an axisymmetric contraction defined by l lt U alters 

7 A contraction such that all cross sections of the tunnel arq similar, whence Is termed 

axisymmeiTic; the sections need not be circular. 
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these waves by a factor of 1 jl x in the longitudinal velocity com- 
ponent and a factor of Iff in the lateral velocity components. 
These equations apply only to particular types of wave; yet when 
the contraction effect is later integrated over the random aggrega- 
tion of waves representing isotropic turbulence, the over-all 
results are found not to differ greatly from the simple factors 
l/li and l /k, respectively. 

Tlie same factors were obtained by Prandtl (ref. 1) for 
other special disturbances: the factor l/l-i from energy con- 
siderations for purely longitudinal disturbance velocities, and 
tho factor Iff from momentum considerations for a rotating 
cylindrical element of fluid with axis normal to the stream. 

EFFECT OF CONTRACTION ON SPECTRUM AND CORRELATION TENSORS 

Effect of contraction on correlation tensor. — The analysis 
herein leads first to tho changes in the spectrum tensor 
Y a f l (k)-*T a p B (ji)- Then the corresponding changes in the 
correlation tensor may be obtained from the Fourier trans- 
form relation, equation (6b) : 


Rat A (r)=f J J 

03 

(16a) 

CO 

•ffa*> B (£)= J* J* J r^"®e'vIcZT(Y) 

(16b) 


03 


In the Fourier integral representation q a A is to be interpreted 
as dg a A , qff as dq a B , Q A as Q a A dr(k), and Q a B as Q^drif). 
Accordingly, 

< 18 > 


since lilff=dT(k)/dr(jc ), Thus 

Qt> A *Q y A k y k a , Q y A Q, A *k y hk a k e \-\ 

z Y v z Y w )J 

The corresponding relation between the postcontraction 
and precontraction spectrum tensors is, by virtue of equa- 
tions (17), 

r ^-w[ IV ® + 2 

' T y s A (k)k y k s k a k 0 \~\ 

Iflfd )j 

where k is related to k by 


In succeeding paragraphs r^ fl (K) will be determined in terms 
of the initial spectrum tensor T a /(i) for various cases. 

Spectrum tensors at upstream and downstream stations 
in terms of Q a . — In an earlier discussion, the Fourier coeffi- 
cients Q a were chosen so as to define a field of turbulence 
confined to a large parallelepiped of volume r, and vanishing 

8 -jt 3 

everywhere outside; for this case lim Q a Q»* was to be 

r— ♦» T 

identified with tho correlation tensor r a/3 . For station A 
upstream of the contraction it will be convenient to special- 
ize this parallelepiped to a cube of edge D. Such a cube 
will, however, be distorted into a parallelepiped of edges 
Iff, Iff, Iff) by the stream contraction by the time it reaches 
station B downstream. (See fig. 1(b).) The spectrum 
tensors for stations A and B, respectively, are therefore 


kl, k 2 , k 3 — li Ki, liK 2 , l$K 3 (20) 

Special case : axisymmetrie contraction but arbitrary 
initial spectrum. — When the contraction is axisymmetrie 
(h—h), the equation of continuity in the form 

Si ? V =0 (21) 

y 

may be used to simplify equation (19). The result may be 
written 





(J'ai A h+T l ff A kg)ki(l — e) 


Y n A kfk a k B (l-eT\ 
{ekf+kf+kff _ 


( 22 ) 


r o /(£)= lim^ Q a A (®Qe A *(k) 

V (17) 

r „/'(K)=]irn^ r3 Q a B ff)Qff(s)\ 

Evaluation of spectrum tensor at downstream station, 
general case. — The identifications made in the last para- 
graph allow the spectrum tensor to be evaluated at station 
B in terms of the spectrum tensor at station A and the param- 
eters l u l i; l 3 defining the stream contraction between sta- 
tions A and B. For a single plane wave q a A =Qffe < iz, 
which is transformed by the contraction into q a B =Q a B e)-£, 
equations (13) give 



where the ratio e^lf/lf; for a large speed gain in the con- 
traction €<C1. 

Special case: axisymmetrie contraction, isotropic initial 
spectrum. — A further simplification occurs when the tur- 
bulence at station A is isotropic. In that case, the spectrum 
tensor r o p fl (V) downstream of the contraction can be expressed 
explicitly to within an unknown multiplicative factor 6(k). 
This results from the fact (ref. 3) that T a f(k) must then be 
an isotropic second-order tensor; the isotropic property 
together with the continuity condition, equation (21), 
requires IV 1 © to be of the .form 

Tafi A (k) = G(k) (klSap — kjef) (23) 


where 


(1 for a=0 

5a * 10 for 3 
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The right-hand side of equation (22) may be evaluated 
by means of equation (23). The diagonal terms reduce to 
relatively simple forms: 



( tf-kfik * 

PPW+W 


14^=4 6(k) 


t 


k 2 -k 2 2 


2fc 1 2 ly(l — e) 
eki 2 +k 2 2 +ks 2 


(24) 


(ekr+W+k/f _ 


(25) 


aDd r M B (V) is obtained from IV s (V) by replacing k 2 by k 3 and 
vice versa. The ’relation, equation (20), between k and * 
applies here. 

One-dimensional longitudinal spectrums. — If the form of 
the intitial spectrum tensor T^Qc) is known, the correspond- 
ing one-dimensional spectrums F a A (ki) can be calculated, 
according to the defining equation (7) as applied at the 
upstream station A: 



k 2j k 2 )dk 2 dk 3 


(26) 


A particular case of isotropic turbulence is of special interest 
(ref. 4): In equation (23) for r^ A (t) the function G(k) is 
taken to be JVfX^+Y 2 ) -3 , where N,y are constants. Then 


and to define F ^(k^—lc'F a B , such that J F a B (Jki)dk\ — 
J F a B d ki] thus 

03 

F ^-m!f r "(vv%) dk ’ ilc ' (27) 


The spectrum tensor elements V aa B following an axisym- 
metric contraction have been evaluated in equations (24) 
and (25). With these values inserted and G(k) specified as 
before, the integrations of equation (26) are best effected in 
polar coordinates. The results are expressed most simply 
in terms of a ‘normalized’ longitudinal wave number kjy 
as incorporated in the two parameters 


«*f+i 




The final result for the one-dimensional longitudinal spec- 
trum following an axisymmetric contraction (4=4) is 


r , N(h 2 +kf) 

11 (^d-^+^+Y 2 ) 3 ' 


and after integration 


FAh) 


ttN 

XW+y*) 


This one-dimensional longitudinal spectrum is of the same 
form as an empirical relation obtained in reference 5 for that 
of isotropic turbulence in the initial period; this agreement 
is the special virtue of the form assumed for GQc). 

The one-dimensional lateral spectrum functions cor- 
responding to the same 0(k) are readily evaluated; they 
are 


F A {k})=F A {ki) 


T N(Zk t 2 + y 2 ) 
2{k 1 2 +y 3 ) 2 


The equality of the F 2 and the F 3 functions results, of course, 
from the isotropy of the turbulence. 

The effect of the stream contraction on these one-dimen- 
sional spectrums is found by employing the postcontraction 
value of r aa , that is, the value r a a B appropriate to the 
downstream station B. Since r aa B is a function of the local 
wave-number vector k at station B, the equation correspond- 
ing to (26) is 


F a 


-JJ 


r oa fl (*„ K 2 , K 3 )dK 2 dK 3 


which is a function of x t . 

For performing the integration and making later com- 
parisons of spectrums, it is convenient to transform from 
Ki 3 k 2j k 3 to ki 3 k 2 , k 3j where 

Ki=kifli 

K 2 =k 2 jl 2 
' K 3 =k 3 /l3 


(2+4«+2l+»l+ 5i) In (j^j)] (28) 


The corresponding result for the one-dimensional lateral 
spectrums following an axisymmetric contraction is 


W*i)=W*i)= 


irN f (2,8—2)? 


2 T ! 4 ,¥ V s 


( ! 


(1 — e)(«— 1) { 


2s-H i s-\-t 


2s 


t 


'■ In 


Gu) + 


1 — e r6s+5t 

2t L 2 


(29 > 


and for «<C1 (large speed gain) a simple but very close ap- 
proximation is 


W*i) 


ttN l-\-2ki 2 /y 3 


fO(e,«lne) 


(The corresponding approximation for F B (k{) is not simple 
enough to warrant noting.) The parameters 4» 4=4, and 
e=4 2 /4 a are related to the initial and final Mach numbers of 
the main stream, if the fluid is air, by the equations 



>- 


(30a) 



SPECTRUM OF TURBULENCE IN A CONTRACTING STREAM 


107 


For incompressible flow (Mg, M A -*0) these reduce to 

h=U a fU A \ 

k U=h~ m > . (30b) 

e=li~ 3 ) 

These postcontraction spectrums, equations .(28) and 
(29), are compared with the initial spectrums in figures 2 
and 3, respectively; the comparison is based on an assumed 
initial stream Mach number of 0.05 (station A) followed by 
an axisymmetric contraction such that the final Mach num- 
ber is 2.0 (station B); the corresponding parameters are 
Z x =29.8, ^=0.382, €=0.00016. Consider first the longi- 
tudinal spectrums, figure 2. Normalizing factors are in- 
corporated such that the areas under the two curves, if re- 
plotted on a linear scale, would be the same; this normali- 
zation serves to differentiate changes in shape from changes 
in amplitude. The figure exhibits a rather striking distor- 
tion of the spectrum after traversing the stream contraction: 
the peak spectral density is shifted from zero wave number 
to &i/y= 1.4 along with a general shift of density to the 
higher wave numbers. Associated with this change in shape 
is a reduction in amplitude by the factor u a 2 and u B 2 

being the respective integrals of the spectral density curves. 
These integrals are evaluated in a later section. 



A .5 .6 .7.8.9 
Wave number 


| Z 3 4 5 6 7 8 9IO 

parameter, k\/y 


FrntjEB 2. — Comparison ot one-dimensional longitudinal spectrum upstream (i-'f ) and 
downstream ( I'f) of arisymmotrio contraction. Curves normalized to same area. Iso- 
troplo Initial turbulence with O(t) Ma— 0.03, Mb<= 2.0, corresponding to 

li=MAil, «-0.0001637. 


• The corresponding comparison for the lateral one-dimen- 
sional spectrums is made in figure 3. In this case the axi- 
symmetric contraction has made very little distortion in the 
spectrum. There is again a change in magnitude (this time 
an increase) in the ratio v B 2 /v A 2 . 

The changes in magnitude (that is, the changes in area 
under the spectral density curves) correspond to the changes 
ujfju^ and v r rlv A in the mean-square components of turbu- 
lence and are, at least qualitatively, well-known. The pre- 
dicted changes in the shape of the spectrum curves are ap- 
parently new. 

In the above comparisons both precontraction and post- 
contraction spectrums have been expressed in terms of the 
precontraction longitudinal wave number k u whereas the 
local postcontraction wave number is ic x =ki/Zi. Consider, 
however, a representative longitudinal wave which has the 
form cos k t x at station A and cos k x £ at station B. If a: and £ 



Figure 3. — Comparison of one-dimensional lateral spectrum upstream (f£) and down- 
stream (Ff) of axfeymmetrio contraction. Curves normalized to same area. Isotropic 
initial turbulence with M a— 0.05, .Mb” 2.0, corresponding to 29.821 

0.0001 637. 


are identified with the respective distances swept in time t 
by the moving waves over stationary hot-wire probes at 
stations A and B, respectively, then k 1 x=k l U A t and * x £= 



l\ U A t. 


Thus the (temporal) frequency seen by the 


hot wire in both cases is k t U A /2ir. The comparison based 
on k, therefore constitutes, in effect, a comparison of the 
time spectrums that would be seen by stationary hot-wire 
probes, in contradistinction to the space spectrums 
discussed in the earlier parts of the paper. 


EFFECT OF CONTRACTION ON MEAN-SQUARE VELOCITY COMPONENTS 
FOR ISOTROPIC TURBULENCE 

The mean-square velocity components of the turbulent 
field may be identified as the diagonal terms of the corre- 
lation tensor R a p(r) with r set equal to zero. Thus 

v?=R n ( 0 ) 

t?=I?22(0) 


7O 2 =J? 33 (0) 

where u, v, w have been written for g x , q 2 , q 3 , respectively. 
The evaluation of these means is much less laborious than 
the evaluation of the general correlation tensor. In par- 
ticular, the evaluation of the ratio of the means ■ / u s 2 l'n A 2 , 
etc., may be made when the initial turbulence is specified 
to be isotropic' but no further details of its spectrum are 
known. These ratios will be calculated in the 1 following 
paragraphs. 

Evaluations of u ! and a 1 at upstream station. — According 

to equation (6) 

00 

^A0)=^?= J J J I V©dr(£) 


For isotropic turbulence T a $ has the form specified in equa- 
tion (23), whence 

co 

J j j om^-mdrQc) 
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where G(k) is on arbitrary function. It is convenient to 
transform to spherical polar coordinates: 


ki=k cos 0 . 


k 3 =k sin 0 cos <p 
k 3 =k sin 0 sin <p 


y 


(31) 


together with for an axisymmetric contraction, whence 


— i frrw(i)(p-wr© 

* ~h 2 J J J (ek S+kf+kff 


Again the polar-coordinate transformation (equation (31)) 
is made, with the result 


Then 


dr(k)=k 2 sin 0 ddd<pdk J 


u A = J* k*6(k)dk J** d<p sin 3 0d0 


i C" ■ C2x rx 

Ujf—n k*Q(k)dk d<p 
k Jo Jo Jo 


sin 3 OdQ 


(e cos 2 0+sin 2 0) 2 


(32) 


The first two integrals occur also in u A 2 (equation (32)), and 
they cancel in obtaining the ratio u, 2 lu/) thus 


For the present purpose the function G(k), which, together 
with the condition of isotropy, defines the turbulence, may 
be left unspecified; the integral involving G(k) will cancel 
out in for min g the ratio Let this integral have the 

value H] then 

«?=! ^ 


Ub _ 3 f r sin 3 0 dO 
u A 4ld Jo (f cos 2 0+sin 2 0) 2 

The final result may be written 


U B 2 3 r -1 
u? 4d 2 [_1 — « 


— — - — tanli _1 Vl — £ 

(!- e )3/ S 


(34) 


By virtue of the assumed isotropy 


and an asymptotic expansion for small e is 


v A =w A 2 =^ irll 


=== — '777 1+^ e+(l + «)ln ^-+0(6 2 ln e) 
u A 4ti L. 2 4 _ 


Evaluation of ratio of u 1 at downstream station to 1? at 
upstream station. — The mean value u B 2 is obtained from an 
integration involving the spectrum tensor after the latter 
has been transformed by passage of the flow through the 
tunnel contraction; according to equation (16b) 

oo 

B u «( 0)=^7=JJJ iV&dTGc) 

— CO 

For the present case, where the spectrum tensor at station 
A is assumed isotropic and the contraction is axisymmetric, 
the transformed tensor r n B (K) has been determined in 
equation (24) . Thus 

— is r c r 

UB ~h J J J (ekd+kJ+k/f 

Because of the unspecified function G(k), it is convenient to 
change the variables of integration from the components of 
k to the components of k. In other words, a transformation 
is made from the "wave-number space” of station B to the 
“wave-number space” of station A. The transformation 
follows from the Cartesian relations 

dT(k)=dkidkidk z \ 

dr (K)=dKidn 3 dK S \ (33) 

_ dki dk 2 dkA 

“irirx; 


Equation (34) gives the ratio of the mean-square longi- 
tudinal velocity fluctuation downstream of an axisymmotric 
tunnel contraction to the corresponding mean square up- 
stream of the contraction, when the initial turbulence is 
isotropic. The contraction is characterized by an increase 
in the stream speed in the specified ratio li and a decrease 
in the lateral dimensions in the specified ratio Zj; the param- 
eters It, h, and t ^IS/l 2 are completely defined by the initial 
and final Mach numbers of the stream according to equa- 
tions (30a) and (30b). 

The variation of -yJu B 2 /u A 2 with the speed, ratio l\ is plotted 
in figure 4 for two examples: in the first the flow is assumed 
compressible with a Mach number 0.05 at the start of the 
contraction; in the second the flow is assumed incompres- 
sible (M a , M b ^> 0). The Mach number scale at the bottom 
applies only to the compressible case, the k scale to both 
cases. The salient characteristic of the curves is the marked 
reduction in the longitudinal component of turbulence with 
increasing speed ratio 

Compressibility is seen to have but a secondary effect, 
which is appreciable only at supersonic speeds. Note 
(equations (30) apd (34)) that with U os the independent 
variable, the effect of compressibility appears only in tho 
parameter e. The physical significance of e follows from tho 
definition of h as the speed ratio provided by the contraction 
and IS as the area ratio of the contraction (in the axisym- 
metric case considered), with t=lS/h 3 - For supersonic final 
speeds it is more proper to speak, of a converging-diverging 
nozzle than a contraction, the term ‘contraction’ having 
been retained herein primarily for reasons of past usage. 
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The basis of the compressibility effect may be summed up 
in the following way. The influence of an axisymmetric 
stream contraction arises from distortion of the fluid ele- 
ments, as described by the parameters l x and l- 2 . (See fig. 1 (b).) 



Downstream Mach number, M g 


Fioobe 4.— Typical examples of selective effect of axisymmetric stream contraction without 
decay on components of turbulent Intensity, showing influence of compressibility. Iso- 
tropic turbulence at li~l. 

These parameters are related by the continuity condition 
where a- is the density ratio. Thus compressibility, 
in allowing <r to deviate from unity, changes the relation 
between k and l 2 somewhat, and consequently modifies the 
contraction effect. 

The graph of equation (34) in figure 4 is primarily for 
illustrative purposes; a form more useful for engineering 
applications is given in figure 5. The single curve provides 

the variation of -yJu B 2 /uP with both h and «; k and e may 
be determined from the initial and final Mach numbers by 
means of the simple relations (30a). __ 

Evaluation of ratio of V s at downstream station to v* at 
upstream station. — The value of vp results from an inte- 
gration involving the transformed spectrum tensor, accord- 
ing to equation (16b) 

00 

i?« B ( o)=^?= r r fiy^TOc) 


For isotropic initial turbulence and an axisymmetric con- 
traction, the transformed spectrum tensor Tx 2 B (k) has been 
evaluated in equation (25). Thus 


i 


^=*iJJJf?(A-)[* 2 -^ 


2k 1 1 k 3 2 (l — e) , 

ekf+W+kf* 


*,**.*(** — * x *) (i — 


dr(jc) 



Figure 5. — Variation of relative root-mean-squan? longitudinal velocity fluctuation with 
both speed ratio li and area ratio l\ for axisymmetric contraction and isotropic turbulence 
ath™l. 


Again it is convenient to transform from K-space to -space 
(equations (33)) and to introduce polar coordinates k, <p, 
and 0 (equations (31)). The integrations with respect to k 
and <p are readily disposed of, with the result 




sin 3 0 cos 2 Odd 
sin 2 0+e cos 2 0 


, C T sin 5 0 cos 2 6dd ~~| 
e Jo (sin 2 d+e cos 2 0) 2 J 

where j k'0(k)dk, as before. Upon carrying out the 

g 

integration and dividing by v a 2 =^ttH, there is obtained 


finally 


Vb 3 1~2 
v A ‘ _1 


2 — t 


-e (1-e) 3 ' 2 


3 

tanli - 


Vi-.] 


(35) 


For small e this has the asymptotic expansion 


Vb 


ri{[ 


2 + e+| « 2 ln e 


]+0(e 2 )| 


Equation (35) gives the ratio of the mean-square lateral 
velocity fluctuation downstream of an axisymmetric tunnel 
contraction to the corresponding mean square upstream of 
the contraction, where the turbulence has been assumed to 

be isotropic. The variation of -yJVp/v/ with the speed ratio 
h is plotted in figure 4, which already contains the graph of 

-y/ u b 2 /u a 2 discussed earlier; again the two cases are incom- 
pressible flow and compressible flow with an initial Mach 
number of 0.05. For Z L > 1 and incompressible flow, the 
lateral component of turbulence is seen to increase steadily 
with li, in marked contrast to the decrease exhibited by the 
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longitudinal component. The curve (of the lateral com- 
ponent) for compressible flow begins to differ sensibly from 
the curve for incompressible flow for downstream Mach 
numbers above 0.3; above sonic speed compressibility is seen 
to effect a complete reversal of the curve. The over-all 
effect of compressibility on the contraction effect is thus 
much greater for the lateral than for the longitudinal com- 
ponent of the turbulence. 

The graph of equation (35) in figure 4 is primarily illus- 
trative; a form more useful for engineering applications is 
given in figure 6 . The single curve provides the variation 
of -\Ivb 2 /Va with both h and e; l, and « may be determined 
from the initial and final Mach numbers by means of equa- 
tions (30a). 

-DECAY CONSIDERATIONS 


1-5, 



s i i i mini mi i i nimi j i mm] i i milt 

.0001 £>01 .01 .1 ID 10.0 

Contraction parameter, « * l|/l^ 

FmuBE 6— Variation ol relative root-mean -mi uaro lateral volodty fluctuation with both 
speed ratio ii and area ratio l\ for aiisymmotrio contraction and isotropto turbulonco at 
Ii~l. 


CRITERIA FOR NEGLIGIBLY SMALL DECAY 

The basis of the present analysis of the contraction effect 
is embodied in equations ( 11 a) relating the precontraction 
and postcontraction vorticity distributions. The simplicity 
of this result and its derivation arises from the neglect of the 
turbulent decay; by decay is meant the viscous dissipation 
and all the (nonlinear) intermixing processes of the eddies 
which together cause the mean turbulent intensity to dimin- 
ish with time. The postulation of an inviscid fluid elimi- 
nated the viscous dissipation, and the limitation to very 
weak turbulence eliminated the intermixing processes. 
(While there can be no dissipation in an inviscid fluid, the 
intermixing processes ordinarily associated with decay will 
occur.) In order to assess the influence of these assumptions, 
equations ( 11 a) will now be derived in a more general fashion 
with the Navier-Stokes equations as the starting point. For 
simplicity the fluid is taken to be incompressible, since the 
major conclusions are unaffected thereby. 

General formulation of changes in vorticity. — By rearrange- 
ment and cross differentiation to eliminate the pressure term 
(ref. 8 , p. 578), the Navier-Stokes equations can be trans- 
formed into < 


D 


oil 


Dt 


= oi, 


dq/ , d?i'_ 

dxi 


- 0)2 


dx 2 


bft' 

d£T 


- vV 1 ol. 


(36) 


and two similar equations, where o »=b>i, a jj, o>, is the vor- 
ticity and g'=g[ /, ft', q» is the resultant velocity. Now let 
q’ be the sum of a stream velocity U, V, W and a turbulent 
velocity field q_~q,, ft, ft; also, let curl U, V, W== 0 , so that 
a is just curl q. Then equation (36) becomes (in tensor 
notation) 


DW1 ~- a ^ +cc ^ +vV2cCl 


Dt 


(37) 


Contraction Decay 


then the neglected decay terms will be examined and cri- 
teria for their neglect arrived at. 

Neglect of decay terms. — Equation (37) minus the viscous 
terms reads, in expanded form, 

P<* d(Z7+ft) | d(?7+ ft) . d(Z7+ ft) 

+Wa + " 3 ""5i7" (88) 


In this and the earlier equations ^ is the 'Lagrangian’ de- 
rivative following the fluid motion. Now consider a line 
segment 8 x u 8 x 2 , 8 x 3 following the fluid motion: its Lagrangian 
derivative can be shown to be 


D 8 x 1 
Dt 


= 8 x 1 


d(U-\-qi) 

daq 


8x 3 


i>x 2 0x3 


(30) 


and two similar equations. It can be seen that a solution of 
equations (38) and (39), together with then companion 
equations, is given by 

coj, C02, ai3~8xi, 8 x 2, 8X3 

for all time t; this result is well known. Now complete the 
neglect of the decay terms by omitting the terms in ft in 
equation (38) and correspondingly in equation (39). By 
this neglect the turbulent perturbations of the flow stream- 
lines have been suppressed: this can be inferred from the 
revised equation (39). If the particles are at station A at a 
time t = 0 and reach station 3 at time t=t, there results 


o>i B 8 x i B 

8 x{^ 


and two similar equations. 

i 5x 3 b . r 
and ^ m Z 3 . 


But 


8 xi B 

Sx ! A 


is just li, 



is k, 


and there are two similar equations. The first set of terms 
on the right-hand side is identified as the contraction effect, 
the second set as the decay effect. First the decay terms 
will be neglected in an attempt to recover equations ( 11 a); 


Therefore equations (lla) have been recovered for the in- 
compressible case (density ratio <r=l). 

Consideration of inertial decay terms. — In equations (37) 
to (39) the decay terms not involving v are the inertial or 
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intermixing terms. These are seen to be nonlinear, 
condition for their neglect is evidently 


bqi 


bU 



bx? 


The 


(40) 


and two similar conditions between q 3 and V, q s and W, 
respectively. In a contraction like that of a wind tunnel, 

the dominant velocity gradients will be and 

these will be of the same order of absolute magnitude. A 
sufficient condition to replace expression (40) is therefore 


bq a 


bXf, 


-C 


bU 
bx i 


(41) 


that is, all of the turbulent velocity gradients are very much 
less than the axial gradient of the stream velocity. This is 
essentially the assumption underlying the distortion equa- 
tions (10), which led directly to the vorticity changes (11a) 
in Taylor’s method. 

In statistical terms an approximate inference from equa- 
tion (41) is, for isotropic turbulence, 


4, 




bU 


bxi 


But by definition of X this may be written 

-}Jv? 


-< 


bU 


( 42 ; 


Accordingly, viscosity may be neglected for that portion of 
the spectrum which satisfies the inequality 




bu 

bx i 


(43) 


In the “initial period” of decay, if the inertial decay criterion 
(42) is satisfied, the major part of the speetrum will satisfy 
condition (43) . 

HOUGH ESTIMATION OF MUTUAL EFFECTS OF DECAY AND CONTRACTION 

When the decay effects are not negligible compared with 
the contraction effects (see criterion developed in the last 
section) the theoretical basis of the present theory of the 
contraction effect is violated. Because negligible decay is 
more the exception than the rule, there is considerable 
incentive to attempt to apply the theory outside the valid 
range by means of assumptions concerning the simultaneous 
effects of decay and contraction. 

Suppose, now, the decay and the contraction are considered 
to occur alternatively, in small steps, starting from isotropic 
turbulence. Each stream tube is considered to contract 
stepwise: between steps there is decay without contraction; 
at each step there is a sudden contraction without decay. 
Let the change in speed ratio per step be dli, the reduction 
in u 2 due to decay be (dv?) D , and the reduction in u 1 due to 
contraction be (du^a. Express the effect of decay in the 
absence of contraction in the form 


The “microscale” X may be interpreted as a sort of average 
eddy diameter weighted in favor of the smaller eddies. 
Equation (42) may be accepted as a practical criterion for 
the neglect of the inertial decay terms, equivalent to one of 
the two assumptions underlying equation (10). The other 
assumption, neglect of viscosity, is considered next. 

Consideration of viscous decay term. — The viscous decay 
term in equation (37) is the term containing v. This term is 
linear and so will affect individual plane waves separately 
without mutual interference. The magnitude of the term 
may be estimated to a sufficient approximation by considering 
a wave carried along by the contracting stream and neglect- 
ing (for this term only) the distortion of the wave imposed 
by the contraction. Thus a component of the wave may be 
written 


Then, if the inertial decay terms of equation (37) are negli- 


gible, the equation reads, with 
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(4) =Wi) (44) 

\u/J D 

where l x is a function of the time of travel (decay time) f, 
and the effect of contraction in the absence of decay in the 
form 



The corresponding differential forms are 


(du 2 ) D 

u 2 

_£'&) dl 
~ D(h) dLl 

(46) 

(du 3 ) c 
u 2 

_ o %) dl 

1 

(47) 


The assumption is now made that equation (46) applies to 
the decay effect per step and equation (47) to the contraction 
effect per step, the only interaction being in the common 
u 2 . 8 This assumption neglects the tendency of the decay 
process to counteract the anisotropy 'introduced, by the 


du J — 

•It Is known that In the ‘Initial’ period of decay — u*di. Equation (46) amounts 
__ ** »• 

to replacing the — u 5 on the right-hand side by "(u 1 ) dmr only; some defense may be 
made of this approximation, considering the progressive deviation »rom Isotropy. 



112 


REPORT 1113 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


contraction effect. Tlie total effect per step is then 


), mo- 
11 2 La*.) a*i)_ 


dU 


whence upon integration the over-all e'ffect is 

%^C{U)D{k) (48) 

u A 2 ■ 

That is, if the effect of contraction alone is expressed hv 
0(h) (equation (45)) and the effect of decay alone by D(l i) 
(equation (44)), then the joint effect imder the assumption 
is expressed by the product 0(li) D(/i). 

Equation (48) is intended to provide a very rough adjust- 
ment of the theoretical contraction effect 0(h) to account 
for decay. This adjustment will be made in the attempt 
to compare the theory with experimental results in which 
the decay effects are of the same order as the contraction 
effects. 

Equation (48) refers to the longitudinal velocity com- 
ponent u; an equation of the same form is obtained for the 
lateral component v. For both cases the function Dili) is 
taken to be the right-hand side of the empirical decay law 
(ref. 9) 

=D(h) (49) 

Kh) 



In table I are listed, for the three experimental arrange- 
ments, the parameters used in the estimation of the decay 
factor (equations (48) and 49)). In reference 10 the initial 
stream velocity U A and scale of turbulence L A were given. 
In reference 11 the value of U A was given, and the value of 
L a was taken to be 0.05 foot, the only scale mentioned; it 
was not clear, however, whether this value of scale applied 
with or without screens. In reference 12 the value of U A 
was inferred from collateral information and is somewhat 
uncertain; the scale L A was estimated from the dimensions 
of the honeycomb. In all three experimental arrangements 
the initial relative levels of turbulence were specified. The 
decay time t of the turbulence was computed as the time for 
a particle to traverse the contraction; the value arrived at 
for Hall’s data (ref. 12) reflects the uncertainty in the 
assumed U A . 

Root-mean-square longitudinal velocity components. — 
The comparison of the theory, including estimated decay, 
with experiment for the longitudinal component of turbulence 
is given in figure 7. The theoretical curve, in each instance, 
is the product of a value computed for contraction alone, 
v neglecting decay, (obtainable from fig. 5) and a value esti- 
mated for decay alone neglecting contraction. (See equa- 
tions (48) and (49).) The agreement with MacPhail’s data 
and with Hall’s single point can be considered good. The 
agreement with the Dryden-Schubauer point, on the other 
hand, is poor; a slight. improvement would result on correc- 
tion for the spurious contribution of the noise background. 


for isotropic turbulence in the initial period. The decay 
time tih) in the formula is the time required by 'a particle 
of the main stream to pass through the contraction, the 
initial velocity being U A and the final velocity liU A . 

COMPARISON WITH EXPERIMENT 

There appear to have been no experimental investigations 
with which to compare the predicted changes imposed by a 
stream contraction on the spectrum of the turbulence, or 
on the correlation tensor of the turbulence. The available 
experimental data seem to be limited to measurements 
bearing on the changes in the root-mean-square velocity 
components. These data apply, moreover, to conditions 
outside the proper scope of the present theory in that large 
decay effects are present. The experimental data are there- 
fore compared with a crude extension of the theory in which 
the decay is allowed for in first approximation. (See pre- 
ceding section.) 

The most extensive data are those of MacPhail (ref. 10), 
which in effect cover a range of contraction ratios from 
^=1 to Zi=9.65 inasmuch as measurements were made at 
various stations along the contraction. Isolated points for 
particular contraction ratios were obtained from investiga- 
tions made for other purposes by Dryden and Schubauer 
(ref. 11) and by Hall (ref. 12). Only those points were 
chosen for which the initial turbulence was indicated to be 
approximately isotropic. 
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Figure 7. — Comparison of predicted axlsymmetrio contraction effect with experiment for 
longitudinal component of turbulence, with decay allowed for In first approximation. 
Initial Isotropic turbulence assumed. * 
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Root- mean-square lateral velocity components. — Com- 
parison of the theory, again including estimated decay, with 
experiment for the lateral component of turbulence is given 
in figure 8. There is complete disagreement with MacPhail’s 
data and Hall’s single point, and on the other hand, good 
agreement with the Dryden-Schubauer single point. Thus 
there is the curious result that AlacPhail’s and Hall’s data 
agree well with theory for the longitudinal component and 
disagree entirely for the lateral component, whereas the con- 
verse is true for the Dryden-Schubauer data. 

Discussion. — The uncertainty both in the manner of 
estimating the decay effect and in the data (table I) on which 
the estimate was based is still far from sufficient to account 
for the discrepancies between theory and experiment for the 
lateral component of turbulence. The very large amplifica- 
tion found by MacPhail is particularly hard to explain. On 
the other hand, the experimental data of the several ob- 
servers show considerable disagreement, especially when 
differences in decay are allowed for. This disagreement 
would tend to cast doubt on the validity of some of the data; 


T 1 1 1 1 r— ( — 

] I I Experiment 

o MacPhail (ref. 10) 

O Dryden- Schuttauer (ret 1 1) 
□ Hall (ref. 12) 


Theory with estimated decay 

for experiment of MacPhail 

for experiment of Dryden-Schubauer 
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Fioobe 8. — Comparison of predicted axtaymmctric contraction eflect with experiment for 
lateral component of turbulence, with decay allowed tor in first approximation. Initial 
Isotropic turbulence assumed. 


the disagreement may also be in part a consequence, pre- 
dicted by the theory, of possible differences of the initial 
spectrums from each other and from isotropy. 

CONCLUDING REMARKS 

The original aim of this paper was to provide a quantitative 
explanation of the observed changes in in!.' toot-mean -square 
velocity components of the turbulence of a wind-tunnel 
stream after passing through the tunnel contraction. The 
simplifying assumption of negligible decay was made to make 
the analysis tractable, although the decay and contraction 
effects are ordinarily comparable. The analysis on this basis 
disclosed, in addition to the above integrated effects, pro- 
nounced changes in the spectrum of the turbulence. The 
changes in the shape of the spectral density curves, as dis- 
tinguished from over-all changes in amplitude, would appear 
to be considerably less sensitive to modification by decay than 
would the mean-square velocity components. For this 
reason, and because such spectral changes have not pre- 
viously been discussed, the emphasis of the present paper 
has been placed most heavily on these spectral effects. 

In particular, it has been found that the one-dimensional 
longitudinal spectrum for isotropic turbulence exhibits a 
rather interesting change in shape downstream of the con- 
traction; the center of gravity of the curve of spectral density 
versus longitudinal wave number is shifted substantially to 
higher wave numbers, the resulting distortion moving the 
peak of the curve well to the right of its initial position above 
the origin. The distortion is quite pronounced and would 
appear to be readily amenable to experimentalobservation. 

The restrictive assumption of negligible decay largely de- 
feats the original aim of the paper. Nevertheless, for 
practical reasons an attempt has been made to provide a 
crude extension to the theory in which decay is allowed for in 
first approximation. With this approximation the theory 
has been compared with experimental values of the contrac- 
tion effect on the longitudinal and lateral component root- 
mean-square velocity fluctuations. The agreement for the 
longitudinal component is good, whereas there appears to 
be almost complete disagreement for the lateral component , 
the experimental data themselves being in conflict. It is 
perhaps premature to attempt any general conclusion. For 
the present, the theory as augmented by the estimated decay 
effect may be useful in wind-tunnel-design applications. 

It is clear that the tunnel contraction effect on the com- 
ponents of turbulent intensity cannot be represented by fixed 
fractional changes independent of the character of the initial 
turbulence. Instead the separate factors for the longi- 
tudinal and lateral components depend markedly on the 
spectrum of the turbulence. For initial isotropy, however, 
unique factors are predicted that, when decay is neglected, 
are independent of the details of the spectrum. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, August 80, 1951 




114 


REPORT 1113 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


APPENDIX A 


SYMBOLS 


The following notation is used in this report: 

The subscripts 1, 2, 3 refer to a rectangular coordinate 
system with the 1-axis alined with the axis of the main flow 
and directed downstream, the 2-axis directed horizontally, 
and the 3-axis directed vertically. Separate systems are used 
with the origins at stations A and B , respectively. (See fig. 1 .) 
Vector and tensor notations are used interchangeably; for 
example, k=k a =(ky, k 2) k 3 ), where a= 1, 2, or 3 designates 
a vector with components k u k 2l and k 3 . 

C(li) function defined in equation (45) 

D(li ) function defined in equation (44) 

D edge length of cube within which turbulent field 

is defined 

e base of natural logarithms 

F a = F u F 2> or F 3 

Fi one-dimensional longitudinal spectral density 

(see equation (7)) ' 

F 3 ,F 3 one-dimensional lateral spectral densities (see 

equation (7)) 

0(k) function appearing in isotropic spectrum tensor 

H j constant k*G(Jc)dk ^ 

1m imaginary part of 

i= , y=i 

K y , K 3 , K 3 =ky+ ky, k 2 -\-k 2 , & 3 +V, respectively 
k amphtude of k(=^lki‘+k 2 i +k 3 3 ) 

k— k a = (ky, k 2 , k 3 ) wave number vector (station A) 

L lateral scale 

la ly, 1?, (3 

ly stream velocity at station B divided by stream 

velocity at station A (see fig. 1) 

1 2 stream breadth at station B divided by stream 

breadth at station A (see fig. 1) 

1 3 stream height at station B divided by stream 

height at station A (see fig, 1) 

M .Mach number of main stream 

N amphtude of special isotropic spectrum tensor 

(see following equation (26)) 

Q== Q a — (Qy, Qi, Q 3 ) disturbance wave amphtude vec- 

tor 

2= 2<r= (ffi, q_ 2 , gj) disturbance velocity vector 

Radix) correlation tensor (ref. 6; 

Re real part of 

r magnitude of £=Vrp4r?+r? 

r= r a =(r lt r 2 , r 3 ) separation vector of two correlated 

points 

(k 3 \ 

=-A+l ) 


t 

U 

u,v,w= 

X 


x= 
r afi(k) 
y 

e 

e a0r 

e 

K 

K= 


V 

f 


<P 


(28) (= 


.(«* 


•1 )ky* 

Zi 


parameter m equation 
time 

main-stream velocity 
q ly 2 s, 23 disturbance velocity components 
length of wind-tunnel contraction (distance be- 
tween stations A and B) 
used occasionally in place of Xy 
x a = (xy, x 2 , x 3 ) position vector (station A) 
spectrum tensor (ref. 3) 

constant in special isotropic spectrum tensor 
(see following equation (26)) (y=l/£) 
contraction parameter (=l 2 /ly i ; see fig. 1) 
alternating tensor defined after equation (9) 
polar angle (equation (31)) 
magnitude of k(=V'c?'+Mi 5 +k3 5 ) 

*a— (m, * 2 , * 3 ) transformed wave number vector 
(station B) ; ( K a =kjla ) 
kinematic viscosity 
magnitude of | ( = VSi 2 + f + fs 2 ) 

£«=(£i» £2, £3) transformed position vector 
(station B) (see equation (10) and fig. 1(b)) 
summation over a for a=l, 2, 3 

stream density at station B divided by stream 
density at station A 
a volume 

azimuth angle (equation (31)) 

G«= (Gi, Gj, G,) vector amphtude of vorticity wave 
io,= (ojy, a^, 0 ) 3 ) vorticity vector 


Superscripts: 

A measured in vicinity of station A, upstream of 

contraction 

B measured in vicinity of station B, downstream 

of contraction 

* complex conjugate 


Subscripts: 
A 

B 

“Ar, 5 

1,2,3 


measured in vicinity of station A, upstream of 
contraction 

measured in vicinity of station B, downstream 
of contraction 

take on values 1, 2, or 3 ’and designate tensor 
quantities 

specific values of a, 0, y, or 6 


A symbol with the mark ~ above it refers to a single plane 
wave. A bar over a symbol designates an average (usually 
a spatial average) ; a bar under a symbol designates a vector. 
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APPENDIX B 


EQUIVALENCE OF SPACE AND TIME AVERAGES IN STATISTICALLY STEADY, HOMOGENEOUS TURBULENCE 


The definitions of statistical homogeneity and statistical 
time independence will first be made precise. Let F (x, y, z, t) 
be some property of a turbulent field that varies ip time 
and from point to point; thus F may be the pressure, or any 
of the velocity components, or a correlation of velocity com- 
ponents at two points of fixed separation, ( x , y, z) being one 
of the two points. If, for all choices of the property F, (a) 
the average of F over a time T—> oo is independent of (x, y, z), 
the turbulence is defined to be statistically homogeneous; 
if (b) the average of F over a volume V— » a> is independent 
of t, the turbulence is defined, in the sense used herein, to be 
statistically steady or time-independent. The respective 
averages are supposed to be approached uniformly, in the 
mathematical sense, as T or V, respectively, approaches 
infinity. (A statistically steady or “stationary” condition is 
defined differently in the theory of random processes.) 

It will now be proved that if the turbulence satisfies the 
two conditions (a) and (b), the time and space averages 
defined therein are equal. In this proof, no resort will be 
made to the “ergodic hypothesis” of statistical mechanics, 
which leads to the equivalence of the time average and the 
“ensemble” average. The possibility of the joint existence of 
the conditions (a) and (b) probably amounts, however, to 
just as fundamental an assumption. 

The space average will be mode over a parallelepiped of 
edges a, b, c and the time average over a time T, and then a 
limiting process will be applied. The average of F over both 
space and time is thus 


_ 1 fTrcrbr'a 

F,',= lim — r-rp Fdxdydzdt (Bl) 

dOCJ. Jo Jo Jo Jo 

Any order of integration is permissible, since the integra- 
tion limits are constants. If the time integration is per- 
formed first the expression may be written 


F,.,= lim —rzrn 

a,b,c,T-><*> CLOCJ. 


rc pb pa / pi 

Jo Jo Jo \J 0 


Fdt) dxdydz 


/ 


By virtue of the postulated uniform convergence of the time 
and space averages, the operation J™may be brought under 
the integral sign: 


F.,,= lim - 1 - f f f riim ^ f Fdt\ dxdydz 
a,b,o-*^o aOC Jo JO Jo \T-m» JO / 

h m: T4xdviz 


-- lim — 

aj),c-+ co CVuC 


(B2) 


where F t is the time average of F. But, by condition (a), 
F t is independent of x, y, and z. Therefore 

(B3) 


Alternatively, the space integration and limiting process 
may be performed first: 


F.,,= lim ~ \ T ( lim ~ f C P f° F dxdydz ) 

T— ►<» -L Jo \a,b,e — *<® (LOC Jo Jo J 0 / 

= lim i f F,dt 

T->cc I Jo 


dt 


(B4) 


where F, is the space average of F. By condition (b), ~F, is 
independent of t; therefore 

F,' t =~F, (B5) 


Equations (B3) and (B5) together state that 

7,=F t =F,, t (B6) 

or the space average, the time average, and the space-dime 
average are all equal. 
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TABLE 1.— DATA FOR ESTIMATION OF DECAY 


Source 

h 

Ujl, 

ft/sec 

L if 

U 

sec 

Vtut 

ft/SCO 

MacPhail (ref. 10) _ 

1.20 

3.55 

0.012 

0.19 

0. 149 


1.60 

3.55 

.012 

.34 



2.55 

3.55 

.012 

.45 



4.90 

3.55 

.012 

.51 



9.65 

3.55 

.012 

.57 


Dryden-Schubauer (ref. 11) 

6.0 

6.86 

.05 

1.31 

.114 

Hall (ref. 12) 

5.2 

1.54 

.025 

1.22 

.048 


F.. t =F, 





